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ABSTRACT 

Stellar jets are normally constituted by chains of knots with some periodicity in their spatial distri- 
bution, corresponding to a variability of order of several years in the ejection from the protostar/disk 
system. A widely accepted theory for the presence of knots is related to the generation of internal 
working surfaces due to variations in the jet ejection velocity. In this paper we study the effect of 
variations in the inner disk-wind radius on the jet ejection velocity. We show that a small variation 
in the inner disk-wind radius produce a variation in the jet velocity large enough to generate the 
observed knots. We also show that the variation in the inner radius may be related to a variation of 
the stellar magnetic field. 

Subject headings: accretion, accretion disks - ISM: Herbig-Haro objects - ISM: jets and outflows - 
stars: magnetic fields - stars; pre-main sequence - winds, outflows. 



1. INTRODUCTION 

Stellar jets are observed in the form of chains of emit- 
ting nebulae called Herbig-Haro (HH) objects. HH ob- 
jects are generally thought to be shock-heated density 
condensations traveling along the outflows formed by 
star- disk systems during th e process of star formation 
(e.g. iReinurth fc Bali^[2nnih . 

Since their discovery several scenarios have been sug- 
gested for the origin of steady outflows from young stellar 
objects. In the stellar wind model, m aterial is acc elerated 
by thermal pressure gradients (e.g. ICantolirQSOf ). Mag- 
netohydrodynamics (MHD) models rely on the magne- 
tocentrifugal launching mechanism ([Blandford &: Pavne I 
[1982). For the X-wind scenario the jet is magnetically 
driven from the so-called "X-annulus" where the young 
star's magnetosphere interacts with the disk (|Shu et al.l 
|2000( ). In the disk wind scenario the jet is launched 
from an extended region of the disk surface (e.g. lFerreira I 
[1997). The analytical models mentioned above have fo- 
cused mainly on the steady-state aspect of the ejection 
phenomena. 

Unsteady periodic ejections with timcscales of the 
order of several rotation periods of the inner disk 
radi us have been obtained by numerical simulations 
(e.e. lOuved fc Pudritd [l997l : iCoodson fc Wingleel [l999l : 
iMatt et al.ll2002f ). Observations indicate a significantly 
longer timescale is associated with the appearance of 
knots in stellar outflows. Nearly all observed jets 
present small scale knots up to 0.1 pc from the cen- 
tral source, with a spacing between the knots corre- 
sponding to a time scale of ^ 1-20 yr (e.g. HH 30 - 
iBurrows et all[l99fit HHlll -iReinurth et al.l [l99l RW- 
Aur - iLopez-Martm et al1l2003l ). More fragmented knots 
are observed on a typical timescale of ~ 10^""^ yr at larger 
distances from the source. While the long term variation 
of jets can be explained by variations in the accretion 
rates (e.g. during FU-Orionis phases), the possible ori- 
gin of small scale knots is still unclear. 



A first possibility, suggested by similarities with extra- 
galactic jets, is that the knots may be fori ned by hydro- 
dyna mics Kelvin-Helmholtz instabilities (jMicono et al.l 
Il998f ). MH D Kelvin-Helmholtz reflective pinch mod e in- 
stabilities (jCerqueira &: de Gotiveia Dal P ino' 1999"). or 
by current-driven instabilities ([Frank et a l. 2000). Nu- 
merical simulations by these authors have shown that 
the shocks generated by plasma instabilities are weaker 
than those seen in observations once radiative cooling is 
taken into account. M oreover, some jets (e.g. HH212 - 
IZinnecker et al.|[l998f ) show a remarkable symmetry on 
both sides of the central star-disk system. This may 
be difficult to explain assuming that the instabilities are 
triggered by small scale perturbations. Additionally, and 
more importantly, optical images in many cases show 
that the compact knots have a bow-shock structure (e.g. 
HHlll - iReipurth et al.lll992D that is difficult to obtain 
supposing that the knots are formed by instabilities. 

A widely accepted theory for the presence of knots in 
stellar jets is related to the generation of internal work- 
ing surfaces due to supersonic v ariations of the ejection 
velocity at the base of the jet ([Raga et al.l [1990). In- 
deed, numerical simulations, using a velocity variations 
of about 10-20% of the average velocity, have been able 
to reproduce the morphology and the emission property 
of the knots in HH objects (e.g. Esquivel et al. 2007). 

In this paper we show that velocity variations which 
lead to the creation of knots similar to the observed ones 
may be generated by variations of the inner disk-wind 
radius. Furthermore, we suggest that these variations 
may be related to changes in the stellar magnetic field. 

This paper is organized as follows. In Section [2| we 
determine the effect of a variation in the inner disk radius 
on the jet velocity. In Section [3| we suggest that the 
variation in the inner radius could be connected to a 
change of the stellar magnetic field. In Section [4| we 
discuss the approximation used and the limitations of 
our model. Conclusions are given in Section [5| 
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2. THE EFFECT OF A VARIABLE INNER EJECTION 
RADIUS ON THE JET VELOCITY 

We first discuss some asymptotic properties of a single 
fieldline ancliored to an accretion disk as a function of the 
conserved invariants and the fieldhne's footpoint radius. 
Then we extend the discussion to an ensemble of field- 
lines anchored at a range of footpoint radii and derive 
mean values for the asymptotic velocity of a disk-wind 
extending from an inner to an outer finite radius. Lastly, 
we discuss the consequences of a time-dependent inner 
disk-wind radius on the average velocity of the outflow. 

2.1. Asymptotics along a single fieldline 

In steady state, axisymmetric, ideal MHD a number of 
quantities are conserved along a fieldline, i. e. a surface 
of constant magnetic flux ^. Among these invariants 
are the mass-to-magnetic flux ratio fc(5'), the total an- 
gular momentum and the corotation frequency or 
ang ular velocity r^( ^) (e .g. JBlandford & Payne 1982; 



Pelletier fc Pudritd 119921 ). These invariants are defined 



as 



~ fc(*) 
fc(*) 



(1) 



(2) 



(3) 



where p is the mass density, Vp, V^, i?p, are the 
poloidal and toroidal components of the velocity and 
magnetic field, and R is the distance from the star. The 
invariants L{^) and the footpoint Ro{^) 

are assumed to be known or given functions. 

L('I') includes contributions from the twisted magnetic 
field. The relation 

L{-^)^Rl{^)n{-^) (4) 

holds, where Ra{^) is the Alfven radius. Eq. ([4]) deter- 
mines Ra{^) and subsequently A(5') = Ra{^)/Ro{^) 
(the magnetic lever arm) completely. Here, and in the 
following, all quantities with an explicitly indicated func- 
tional dependency on 4*, e.g. i(^'), are taken to depend 
on the magnetic flux \1/ only. Similarly, all quantities 
with a subscript o are taken to be defined at the base of 
the fieldline, i. e. the equatorial plane. The corotation 
frequency ^{^) can be easily evaluated at the equator - 
where the poloidal velocity component Vp vanishes - and 
equals the orbital frequency at the equatorial plane f2o: 

r!(*)=0o(*). (5) 
Another conserved quantity is given by the total en- 
ergy E{^) per unit mass, with contributions from ki- 
netic, thermal, gravitational energy $ and the energy of 
the electromagnetic field (Poynting flux): 



^ ' 2 
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P 



n{^)RB^ 



7 — 1 p k(w) 

where P is the thermal pressure and 7 is the ratio of spe- 
cific heats. The asymptotic jet velocity along a fieldline 
5* anchored at the footpoint radius i?o(^) can be esti- 
mated from the Bernoulli equation for large distances 
and negligible enthalpy as (|Michel II1969D : 



,{^) - %/217(^')i?o(*)A(*), 



(7) 
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Fig. 1. — Average outflow velocity as a function of Rinl R\n,o- 
The average outflow velocity is normalized to (vq) = 

The asymptotic jet velocity varies from fieldline' to 
fieldline. However, we are ultimately interested in a typ- 
ical value representative of the jet as a whole. Therefore, 
we will weigh the asymptotic velocity with the mass-flux 
carried by the fieldline. 

The asymptotic mass-flux along a fieldline per unit 
magnetic fiux (for one hemisphere) can be estimated by 
rewriting eq. ([T]) as 



dM 



k{^) 



(8) 



where M is the wind mass-loss rate. 



2.2. 



Asymptotic mean jet velocity 

The jet is assumed to originate in a disk-wind. Mass 
loading onto the disk- wind shall be efficient only for mag- 
netic flux surfaces in the range \l/in < ^' < Vfex, where the 
flux surfaces ^"111 and 4'ex are anchored at i?in = i?o(^in) 
and i?ox = ^o(^cx)j respectively. The mean or typi- 
cal velocity of the jet given by the superposition of the 
asymptotic velocities along the fieldlines anchored be- 
tween i?i„ and i?ex can be determined weighting the jet 
velocity by the mass-flux carried along every fleldline: 



4M 



(9) 



So far, we have assumed steady state and axisym- 
metry. In the following we will add two assumptions, 
namely keplerian rotation and self-similarity (see e.g. 
IVlahakis fc Tsingano£l ll998V 

Assuming keplerian rotation for the plasma in the 
equatorial plane fixes the corotation frequency to 



C!(*) = 17o(i?o(*)) 



/ GM 



(10) 



The magnetic fiux is a power-law infi, i. e. ^'^fJ^". 
Self-similarity fixes the po wer-law index to a = 1/2 
(|Blandford fc Pavne I Il982f ) . In terms of the footpoint 
radius i?0 5 ^ is given by 



«' = *cx {Ro/Ro.?^*. 



(11) 



Self-similarity also fixes the mass-to-magnetic fiux 
function A;(^') to 

;x(i?0/-Rcx) ^^"^ , 



fc(*) 



(12) 
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Fig. 2. — Schematic representation of tiie studied model, where 
an extended disk-wind is ejected by an accretion disk with a time- 
varying stellar magnetic field. 

and renders the lever arm independent of the magnetic 
flux surface (i. e. A(^') = A). Keplerian rotation and con- 
stant A is sufficient to completely determine the asymp- 
totic velocity along a fieldline as 



I 2GM 



(13) 



while the wind mass- loss rate is determined by eq. |5]and 
yields 

dM = ^fccx*ex^ . (14) 

Finally, defining x — Rc^/Rin, the average velocity (eq. 
O can be integrated, giving: 



^1/2 - 1 

(v) ^ 2 Voo(-Rox) • 

Inx 



(15) 



where i?ex is assumed fixed. The velocity variation as 
function of the inner radius, for different values of X: 
is shown in Fig. [TJ As this Figure clearly shows, an 
increase in the inner radius causes a decrease in the av- 
erage velocity. Moreover, the velocity variation is nearly 
independent of the ratio xo = Rex/Rin,o (where i?in,o is 
a position of the inner radius chosen as a reference). To 
produce a variation of 20% in velocity it is necessary to 
have a variation of order of 30% in the inner radius. Fig. 
[T] shows also the difference between the velocity of the 
disk-wind at the inner radius (keplerian, see eq. I13p and 
the average velocity (defined by eq. [Q]). 

3. THE EFFECT OF A MAGNETIC FIELD CYCLE ON THE 
INNER DISK- WIND RADIUS 

So far, we have dealt exclusively with the magnetic field 
of the disk-wind. These magnetic flux surfaces carrying 
the disk-wind are necessarily open to infinity. A non- 
open fieldline cannot contribute to the global outflow. In 
the following we consider the superposition of a dipolar 
magnetosphere carried by the central star, contributing 
magnetic flux \1/dp and the previously discussed magnetic 
held ^E" which threads the disk and carries the outflow. 
Then the total magnetic flux is given hy + 

This conflguration leads to three different kinds of 
magnetic field lines as is illustrated in Fig. [51 Firstly, 
fieldlines anchored in the polar region of the stellar sur- 
face. While these fieldlines are open to infinity, they are 
assumed to carry only negligible mass-flux and do not 
contribute significantly to the outflow. Secondly, closed 
dipolar fieldlines anchored at lower latitudes of the stel- 
lar surface do not contribute to the outflow in our model. 
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Fig. 3. — Variation of the velocity (normalized to the average 
velocity) as a function of 5 (the amplitude of the stellar magnetic 
field variation) for different values of xo- 

Lastly, a global open magnetic field threading the disk, 
which carries the disk- wind. 

The transition from the closed dipolar fieldlines to the 
open disk-wind, i. e. the location of the innermost open 
flux surface ^^in in the equatorial plane, is given by the 
location of the saddle point of the magnetic flux distri- 
bution 



* 



*DP) 



0, 



(16) 



where the dipolar field on the equatorial plane is given 
by a time-varying dipole moment m(<) as 



(17) 



and ^E" is given by eq. [TTJ Since both magnetic flux 
distributions are known, the inner radius of the disk-wind 
is given in terms of the time- varying stellar magnetic field 
(defined on the stellar surface) from eq. (TH] as 



(18) 



For simplicity, the stellar magnetic dipole moment is 
assumed to vary in time as 



m(t) ~ mo {1 + S sm{2TTt / T^,)) 



(19) 



where and S are the period and the amplitude of the 
stellar magnetic field variation, and ttiq is the stellar mag- 
netic dipole moment at t = 0. 

The inner radius of the ejecting region depends on 
the variation in the magnetic field intensity, and changes 
on timescales of t^,. The average velocity will therefore 
change periodically following the magnetic field cyclic 
variation. 

A plot of the velocity variations, that is the difference 
between the maximum and minimum velocity, and nor- 
malized to the average velocity during one cycle of the 
stellar magnetic field, is shown in Fig. [3] as function of 
S. The different curves are calculated corresponding to 
different values of xo (= Rcx/Rin,o)- A 50% variation in 
the stellar magnetic field produces a ~ 20% jet velocity 
variation, nearly independent of the value of xo- 

We have assumed that the only effect of a stellar mag- 
netic field variation on the ejected jet is a change in i?in. 
Actually, some of the open stellar flux will thread the 
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disk. Therefore, the disk magnetic field (including also 
the contribution from the open stellar magnetic flux) will 

not decrease as ~ i?Q and will also change with time 
following the change in the stellar dipole. 

A change in the stellar magnetic field produces changes 
in the mass-flux and velocity. In fact, the lever arm A 
is a function of the magnetic field and of the mass fiux 
(^at is also a fun ction of the magnetic field). iFerreiral 
11997 ) and Cass e fc Ferreiral (|2000,) studied models of 
self-similar MHD accretion disks driving jets. These au- 
thors showed that ~ 1-1-1/ (2^), where ^ is the ejection 
parameter (that relates the mass flux to the radius by 
M ~ r^). Additionally, ^ ~ 0.1^'^, where /i is the mag- 
netization parameter (given by ^ = B"^ /{%np)^ where p 
is the thermal pressure). 

In this case, eq. [T3l and [Ml mav be written as: 



v^{-^) = A(*) 



^MocCi^g— (20) 



An increase of the magnetic field in the disk (due to the 
increased stellar contribution) produces an increase in 
the ejection parameter ^ and of the mass flux, and a drop 
in the asymptotic velocity. Therefore, the jet velocity 
decreases further with respect to the value calculated by 
eq. m 

The average velocity (eq. ^ is given by: 



jl,^ t;oo(i?cx)A/A(i?ex)?X«-3/2da: 



(21) 



where in the integrals x — Rq/Rc^- Assuming ^ indepen- 
dent of X (i. e. self-similarity) , eq. [21] leads to: 



(v) 



1/2 



1-2C x« 



1 



Woo (^o 



(22) 



where R 

cx is cissumed to be fixed-, but Vcxd {Rcx) varies 
with the stehar magnetic field as a function of A(i?cx) 
(~ 1/(20^/2 for C«l). 

For ^ ^ (assuming that A is converging to a large 
but finite value) eq. [22] reduces to eq. [T5l 

The ratio between the external and inner radius x is 
difficult to determine theoretically, but has been recently 
constrained by observations. In fact, interpreting the ob- 
served transverse velocity shifts in T Tauri microjets as 
indication of rotation, a range of ejecti ng radii corre- 
spond ing to X ~ 10 can be inferred (e.g. IFerreira et al.l 
[20061) . 

Therefore, we may assume ^ <C 1, and write the av- 
erage outfiow velocity normalized to (vq) = {v)fi-R^^g 
as: 



(«o> 



(23) 



where the scaling (v) ~ l/-\/C comes directly from Vao ~ 
A ^ 1/ v^. From this equation, it becomes clear that 
the stellar magnetic field variation necessary to produce 
a 20% variation in the average velocity may be much 
smaller than the 50% value determined from Fig. [31 The 
50% value represents an upper limit to the necessary stel- 
lar magnetic field variation. 



4. DISCUSSION 

Is this section we discuss in some detail the approxi- 
mations and the limitations of the model. First, to use 
the results from the standard disk-wind theory, we ap- 
proximated the evolution of the system as a sequence of 
steady-state configurations. 

We made two key assumptions on the behavior of the 
stellar magnetic field. The first assumption of a dipo- 
lar stellar magnetic field at the inner edge of the disk 
is standard in the magne tospheric accretion model for 
young stars (|Konig]||199lD . On the stellar surface, mag- 
netic fields as larg e as several kG have been observed (e.g. 
lJohns-Krul]||1997l ). The dipole component (- 100 - 500 
G on the stellar surface) , which decays more slowly then 
the multipole components, dominates at large distances 
from the star. 

The second hypothesis is that a dynamo mecha- 
nism operates in protostars ([Chabrier fc Kiike"i1 120061 : 
iDobler et al.ll2006l ) with a periodic or quasi-periodic vari- 
ation of the magnetic field with a typical timescale of a 
few to tens of years. A dynamo mechanism is neces- 
sary to explain the observed magnetic field, because this 
would be dissipated on scales of order of i?*//; « 100 
yr, where i?^ is the stell ar radius and rj is the tu rbulent 
magnetic diffusivity (e.g. lChabrier fc Kiikerl[2006[ ) . A dy- 
namo mechanism, if present, is not neces sarily cyclic as 
neede d by our model. However, recentlv [Sokoloff et al.l 
(|2008f l proposed models of dynamos in low mass, fully 
convective stars with a cyclic magnetic field. The 
photometric variability observed in T Tauri st ars (e.g. 
iMel'Nikov fc Grankin|[2005tlGrankin et al.ll2007D may be 
relat ed to a correspon ding change in the stellar magnetic 
field (lArmitagel[l995l ). 

Finally, our results depend on the complex interaction 
between the stellar magnetic field and the disk. This has 
been studied by a numbe r of authors bo th analytically 
and numerically (see, e.g.. IUzdenskvl[2b04l and references 
therein). It is well understood that the differential rota- 
tion in the disk produces a winding-up of the magnetic 
field lines on a timescale r = 2tt/ {fl^ — fl) (where il^ 
is the stellar rotation frequency) with a following open- 
ing of the stellar magnetic field lines and a consequent 
outfiow. 

The later evolution of the system is unclear. Once 
the magnetic field lines are o pened, they can stay in 
that configuration indefinitely (|Lovelace et al.|[l995( ) al- 
lowing a disk- wind to be ejected from an inner radius 



i?in to an outer radius Re 



In this case the stellar 



magnetic field lines are connected to the disk only in 
a small region around i?in (see Fig. [5]). i?in can be deter- 
mined as the distance from the star where the magnetic 
stress B^BpR^ begin to dominates over the viscous stress 

Md{p.Rl) /dRo fe.g. [Wa^[l996l ): 



Ri, 



2Rl 



-1 2/7 



,1/2 



(24) 



where i?* and are the star radius and mass, = 
(GM/R^y^"^, and we assume B^ « Bp (with Bp (x m, 
defined by eq. [T9|l and a full penetration of the vertical 
magnetic field lines into the disk. The same dependence 
of i?in on Bp is recovered and the results obtained in the 
previous section are unchanged. 
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Alternatively the stellar-disk magnetic field lines 
might close again via magne tic reconnection producing 
quasi-periodic eje ctions (e.g. iGoodson fc Wingleel[l999t 
iMatt et all |2002[ ). Our results correspond to the hy- 
pothesis that the mass flux carried by the quasi-periodic 
ejections is much less then the mass flux carried by the 
disk-wind. If on the other side the mass flux carried by 
the coronal mass ejection is important, the stellar cycle 
would also produce a variation in the jet velocity. Ac- 
tually, in contrast with the disk-wind case, for a coronal 
mass ejection an increase in the magnetic field would 
produce an in crease in the eject ion velocity (as found in 
simulations bv lMatt et al.l[200^ . In fact, in coronal mass 
ejections the magnetic energy is directly transformed in 
kinetic energy by reconnection events. 

We concentrated on a cyclic magnetic field as origin of 
the variations in the inner launching radius for the disk 
wind. However, any physical process that leads to quasi- 
periodic variations of the launching region will result in 
similar variations of the mean asymptotic velocity - and 
therefore to the appearance of knots - as long as those 
variations are large enough. 

It is unclear which process (other than a stellar mag- 
netic field vari ation) might create var iations on a '--^ 1 — 20 
yr timescale. ISpruit fc TaamI (|1993D showed that a re- 
laxation oscillator may produce a change in the inner 
radius in X-ray binaries. The timescale of the oscilla- 
tion depends on the details of the stellar magnetosphere 
and can vary by two orders of magnitude. A variation in 
the accretion rate may also produce a similar variation 
in the inner radius ( see eq. and the discussion by 
iBlackman et al.ll200lD. 

iHartigan et al. showed that short period T 

Tauri binaries do not seem to have jet properties differ- 
ent from those of wider binaries or single stars, arguing 
in this way aga inst b inarity as the origin of jet knots. 
IHartigan et al.l (j2004f ) also argued that there is no cor- 
relation between the ejection of a new knot in CW Tau 
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and an increase in source brightness. One might have 
expected this if the knot was caused by an accretion out- 
burst. Finally, iFerreira et al.l (|2006f ) suggested that the 
alternation between X- and Y-type interactions between 
stellar and disk magnetic fields as due to a magnetic field 
cycle could explain the observed periodicity. 

5. CONCLUSIONS 

In this paper we explored the effect of a stellar mag- 
netic field variation on the ejection velocity of proto- 
stellar jets. While large scale knots may successfully 
be explained by a large increase in the accretion rate, 
we showed that small scale knots may be produced by 
changes in the inner radius. 

We showed that a stellar magnetic field variation, if 
present, represents a natural candidate to produce stellar 
knots similar to the observed. In fact, a stellar magnetic 
field variation produces a change in the inner radius and 
therefore in the jet velocity. We also estimated that a 
stellar magnetic field variation < 50% produces a vari- 
ation (^ 20%) in the average velocity large enough to 
produce knots similar to the observed. 

The stellar magnetic field periodic or quasi-periodic 
variation remains an hypothesis of our scenario, but fu- 
ture large scale temporal observations of magnetic field 
in young stars, and progress in dynamo theory, will both 
help to provide an answer to this problem. 
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